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STATEMENT  OF  PROBLEM  STUDIED  (OVERVIEW) 


Network  polymers  constitute  about  25%  of  all  polymer  products,  the  proportion  in  high 
performance  applications  -  high  strength  composites,  protective  coatings,  electronic  encapsulants, 
medical  implants  -  being  even  higher.  But  because  they  are  insoluble  and  their  molecular  structure 
is  more  compl  Jian  linear  polymers,  molecular  characterization  of  network  polymers  is  much 
more  difficult.  The  understanding  of  relations  between  network  polymer  structure  and  properties 
is  correspondingly  less  certain. 

The  first  step  in  understanding  network  polymers  lay  in  the  analysis  of  network  formed 
under  ideal  conditions:  stepwise  systems  with  equal  reactivity,  no  substitution  effect,  no  secondary 
reactions,  no  loop  formation,  and  no  vitrification  or  phase  separation.  In  our  work,  model  silicone 
and  urethane  systems  were  chosen  as  most  nearly  attaining  these  ideals.  The  corresponding 
theoretical  analysis  (branching  theory)  was  based  on  the  recursive  nature  of  the  network  polymers, 
which  permits  relatively  simple  derivation  of  average  structural  parameters  such  as  crosslink 
density,  sol  fraction,  dangling  chains,  gel  point,  and  molecular  weight  before  gelation.  Our 
success  with  ideal  networks,  both  experimentally  and  theoretically,  inspired  us  to  attack  the  more 
complex  problems  of  networks  formed  under  nonideal  conditions. 

One  important  class  of  networks  are  those  formed  from  epoxides  reacting  with  diamines, 
which  are  used  extensively  in  aerospace  applications  and  in  electronics.  These  systems  may 
exhibit  several  nonidealities,  such  as  substitution  effects  (among  the  amino-hydrogens)  and 
secondary  reactions  (etherification).  Another  complication  is  diffusion  control  of  the  reactions  as 
the  glass  transition  temperature  exceeds  the  cure  temperature.  These  problems  translate  into 
difficulties  in  determining  the  structure  of  the  network.  In  the  case  of  substitution  effects  and 
secondary  reactions,  correct  determination  of  the  structure  entails  characterization  of  the  reactions 
combined  with  branching  theory  appropriately  modified  to  include  those  nonidealities.  In  the  case 
of  vitrification,  an  understanding  of  how  the  glass  transition  temperature  changes  with  conversion 
is  needed  before  the  appropriate  branching  theory  can  be  employed.  Once  the  structure  of  the 
network  is  defined,  structure-property  relations  can  be  used  or  correlations  developed. 

Another  important  class  of  networks  are  those  formed  by  free  radical  polymerization,  a  great 
fraction  of  those  being  used  in  molding  applications,  composites,  radiation  cured  coatings,  and  so 
forth.  These  systems  present  a  host  of  nonidealities:  the  number  of  reactions  (initiation, 
propagation,  termination  (both  coupling  and  disproportionation),  and  transfer),  nonideal 
copolymerization,  vitrification,  phase  separation,  Trommsdorff  effect,  and  formation  of  cycles.  If 
one  looks  first  at  cyclization,  the  characterization  of  the  reactions  are  again  the  important  first  step, 
and  it  is  complicated  by  these  intramolecular  reactions  which  are  difficult  to  analyze  directly.  Once 
this  characterization  is  performed,  the  recursive  analysis  can  be  employed  to  determine  the 
structure,  which  again  can  be  related  to  property  via  relationships  or  correlations. 

The  above  is  certainly  not  an  exhaustive  list  of  common  nonidealities.  Others  worthy  of 
study  include  the  crosslinking  of  chains  of  arbitrary  length  and  site  distribution  (in  vulcanization, 
for  example),  ring  formation  in  stepwise  systems,  and  networks  formed  by  living  (anionic) 
copolymerizations. 


SUMMARY  OF  RESULTS 


Probability  Models 
Professor  Doug  Miller 

The  following  work  was  accomplished  or  is  in  progress: 

Two  manuscripts  coauthored  by  D.R.  Miller  and  C.W.  Macosko  have  been  tentatively  accepted 
for  publication  by  Journal  of  Polymer  Science.  Phvsics  Edition  subject  to  minor  revisions.  They 
are  entitled  "Molecular  Weight  Relations  for  Crosslinking  of  Chains  with  Length  and  Site 
Distribution"  and  "Network  Parameters  for  Crosslinking  of  Chains  with  Length  and  Site 
Distribution." 

A  model  has  been  developed  for  A3  +  B2  copolymerization  with  cyclization  which  treats  small 
cycles  exactly  (using  kinetic  analysis  assuming  Gaussian  end-to-end  distances)  and  larger  cycles 
approximately  (using  Gordon's  spanning-tree  approximation).  Enrique  Vall6s  is  currently  doing 
some  numerical  calculations  with  it  to  evaluate  it 

A  model  for  anionic  polymerization  to  the  gel  point  has  been  constructed  and  a  program 
written.  No  cyclization  is  permitted.  This  model  is  being  used  by  Lutz  and  co-workers  at 
Strasbourg  to  study  an  aMS  +  DIB  system. 

The  theory  has  been  worked  out  for  a  model  of  anionic  polymerization  (living  polymers 
crosslinking)  which  allows  cyclization.  It  is  based  on  spanning-tree  concepts.  Execution  is  planned 
under  the  new  proposal. 

EPOXY  NETWORKS 

Structure- Property  Relations  for  Model  Epoxy  Networks 
Sue  Ann  Bidstrup,  Ph.D. 

Rheological  and  molecular  changes  during  epoxy-amine  cure  were  monitored.  An 
understanding  of  the  relationship  between  bond  formation,  structure  and  mechanical  properties  has 
been  sought  For  this  study,  two  common  epoxy-amine  systems  have  been  chosen:  diglycidyl 
ether  of  bisphenol  A  (DGEBA)  and  tetraglycidyl  diamino  diphenyl  methane  (TGDDM)  both  cured 
with  diamino  diphenyl  sulfone  (DDS). 

Branching  theory  utilizing  the  recursive  nature  of  the  branching  process  and  elementary  laws  oi 
probability  has  been  developed  to  predict  average  structural  parameters.  The  possibility  of  a 
substitution  effect  among  the  amino-hydrogens  and  the  existence  of  an  etherification  reaction  was 
included  in  the  derivation.  This  model  was  applied  to  the  DGEBA-DDS  and  TGDDM-DDS  system 
to  obtain  predictions  for  a  variety  of  molecular  parameters  (including  weight  average  molecular 
weight  (MU,  weight  average  molecular  weight  of  the  longest  linear  chain  (M^c),  gel  point  and 
sol  fraction)  as  a  function  of  the  extent  of  reaction,  the  stoichiometric  ratio  anathe  ratio  of  rate 
constants  between  the  various  cure  reactions. 


The  predictions  of  the  recursive  models  were  tested  by  stopping  the  reaction  after  a  designated 
time  interval  during  the  cure  and  subsequently  analyzing  for  certain  structural  parameters  and  the 
conversion  of  functional  groups.  Epoxide  conversion  was  monitored  by  differential  scanning 
calorimetry,  Fourier  transform  infra-red  spectroscopy  and  titration.  Weight  average  molecular 
weight  was  measured  during  the  cure  by  static  light  scattering.  The  weight  fraction  of  soluble 
material  at  a  particular  conversion,  gel  point  and  critical  stoichiometric  ratio  were  determined  by 
extraction  experiments.  Excellent  agreement  was  obtained  between  the  structural  analysis  and 
predictions  of  the  recursive  theory  for  the  DGEBA-DDS  systems  with  excess  amine  or  at 
stoichiometric  theory  for  the  DGEBA-DDS  systems  with  excess  amine  or  at  stoichiometric 
concentrations.  The  predictions  of  the  recursive  theory  were  lower  than  the  data  obtained  for  the 


TGDDM-DDS  system.  This  may  be  a  result  of  cyclization  occurring  in  the  system  or  the  limited 
ability  of  the  statistical  model  to  account  for  stochastic  correlations. 

The  viscosity  of  the  reacting  system  was  measured  as  a  function  of  conversion  of  epoxide 
groups.  Weight  average  molecular  weight  (Mw)  and  molecular  weight  of  the  longest  linear  chain 
(Mllc)  predict  from  the  branching  theory  were  considered  as  possible  correlating  parameters  i.  -r 
the  viscosity  of  the  system.  The  friction  factor  C  was  also  calculated  during  the  cure.  Good 
correlation  was  obtained  between  viscosity  and  £mw  for  the  DGEBA-DDS  system  at  balanced 
stoichiometry. 

Glass  Transition  Temperature  in  Edoxv  Networks 
Arturo  Hale 


A  literature  survey  has  shown  that  the  prediction  of  Tg  as  a  function  of  conversion  in  thermo¬ 
setting  systems  is  still  an  open  problem.  The  case  of  epoxy  resins  is  especially  important,  since  in 
several  of  these  systems  the  glass  transition  temperature  may  reach  and  even  pass  the  cure 
temperature,  and  the  reaction  becomes  diffusion  controlled.  Thus,  a  full  description  of  the  kinetics 
of  the  polymerization  reaction  requires  knowledge  of  how  Tg  changes  with  conversion. 

Due  to  the  lack  of  published  data  describing  the  variation  of  T_  with  degree  of  conversion,  we 
have  collected  this  kind  of  data  for  a  DGEBA-DDM  system.  We  nave  also  studied  the  behavior  of 
novolac-epoxies  cured  with  novolacs.  This  system  finds  application  in  the  encapsulation  of 
integrated  circuits.  The  basic  chemistry  has  been  studied  by  DSC  and  NMR  using  model 
compounds.  It  has  been  found  that  no  significant  amount  of  competing  reactions  takes  place  if  the 
concentration  of  catalyst  is  low  enough.  The  dependence  of  T  on  conversion  was  studied  for 
epoxy-novolacs  with  different  functionalities.  T  shows  a  very  steep  increase  at  high  conversions, 
which  could  be  due  to  a  reaction  taking  place  between  secondary  hydroxyls  and  epoxide  groups. 

ACRYLATE  NETWORKS 

Formation-Structure  and  Structure-Property  Relationships  During  Vinvl-Divinvl 
CoDQlvmerization 

Donald  T.  Landin 

The  network  forming  copolymerization  of  vinyl  and  divinyl  monomers  has  been  investigated. 
An  understanding  of  the  relationship  between  bond  formation,  polymer  structure,  and  polymer 
properties  has  been  sought  The  main  focus  of  this  understanding  is  on  the  pre-gel  portion  of  the 
reaction.  The  viscosity  of  the  reaction  mixture  is  the  polymer  property  of  primary  interest  at  this 
time. 

The  particular  chemical  system  studied  consisted  of  methyl  methacrylate  with  small  amounts 
(0-5%  v/v)  of  ethylene  glycol  dimethacrylate.  The  radical  addition  polymerization  was  initiated 
with  azobisisobutyronitrile.  A  transfer  agent  was  used  to  control  the  primary  chain  length. 

Bond  formation  was  studied  experimentally  by  measuring  monomer  conversion  and  pendant 
vinyl  conversion  as  a  function  of  reaction  time.  Monomer  conversion  was  found  to  go  by  first 
order  kinetics  at  low  conversions.  1  H  NMR  analysis  showed  that  in  this  system  the  different 
monomeric  units  go  into  the  polymer  equally.  Pendant  vinyl  conversion  was  measured  by 

NMR.  When  plotted  vs.  monomer  conversion,  it  gives  a  nearly  linear  polot  with  a  positive 
slope  and  a  positive  y-intercept.  The  y-intercept  is  attributed  to  cycles  formed  during  the  formation 
of  a  primary  chain.  The  slope  of  the  plot  is  due  to  crosslinking  and  cycles  formed  after  a  crosslink 
has  formed.  This  concept  of  cyclization  is  incorporated  into  a  kinetic  model  for  the  polymerization 
Rate  constants  are  obtained  from  experimental  data.  Pendant  vinyls  has  a  lower  reactivity  than 
monomeric  vinyls.  The  kinetic  model  can  be  used  to  predict  the  conversion  of  pendant  vinyls  to 
primary  cycles,  secondary  cycles,  and  crosslinking. 


A  branching  theory  utilizing  the  recursive  nature  of  the  branching  process  has  been  developed 
to  predict  polymer  structure.  The  monomer  and  pendant  vinyl  conversions  provide  the  input 
information  for  the  theory.  Expressions  for  weight  average  molecular  weight  (Mw),  weight 
average  molecular  weight  of  the  longest  linear  chain  (M1>w),  weight  fraction  of  solubles,  and 
crosslink  density  are  derived. 

The  viscosity  of  the  reaction  mixture  is  measured  as  a  function  of  time  and  then  related  to 
conversion.  w  correlates  the  viscosity  well  as  the  viscosities  of  different  systems  fall  onto  one 
curve. 


Neil  Dotson 

Theoretical  study  of  departures  from  treelike  structure  in  free-radical,  crosslinking  systems  has 
been  initiated,  concentrating  first  on  cyclization  (intramolecular  reaction  of  pendant  vinyl  groups). 
The  groundwork  for  two  theoretical  approaches,  a  probabilistic  model  and  a  simulation  model,  has 
been  performed. 

The  probabilistic  model  is  based  on  competing  reactions  (ring,  chain)  and  employs  equilibrium 
chain  statistics.  The  work  thus  far  has  been  limited  to  single  cycle  growth  on  an  isolated  primary 
chain,  as  this  should  be  indicative  of  the  extent  of  pendant  conversion  extrapolated  to  zero 
monomer  conversion  at  low  crosslinker  levels,  which  has  been  measured  by  several  researchers, 
including  Landin.  Comparisons  between  theory  and  experiment  show  that  the  model  fairly 
consistently  overpredicts.  Possible  causes  few  this  are  excluded  volume  effects  and  end-to-back- 
bone  effects.  This  model  is  also  incapable  of  predicting  any  crosslinker  concentration  dependence, 
another  weakness  which  must  be  resolved. 

A  simple  Monte  Carlo  simulation  model,  similar  to  that  of  Herrmann,  Stauffer,  and  Landau 
(7.  Phys.  A.,  16, 1221  (1983)),  has  also  been  used  to  study  the  zero-conversion  limit  The 
advantage  of  this  method  is  its  explicit  embeddedness  in  space,  so  that  structures  more  complicated 
than  a  single  cycle  can  be  built  up  naturally.  This  model  has  been  extended  to  include  implicit 
small  molecule  motion,  but  still  lacks  macromolecular  motion.  Future  work  will  address  this 
problem,  as  well  as  the  question  of  the  lattice  constraint  and  of  higher  conversion  situations. 


PUBLICATIONS 


New  (attached  as  appendices) 

Calculation  of  Average  Network  Parameters  Using  Combined  Kinetic  and  Markovian  Analysis 
(D.R.  Miller  and  C.W.  Macosko),  submitted  for  publication  in  "Biological  and  Synthetic 
Networks",  Elsevier,  to  appear  1987. 

"Cyclization  and  Reduced  Reactivity  of  Pendant  Vinyls  During  the  Copolymerization  of  Methyl 
Methacrylate  and  Ethylene  Glycol  Dimethacrylate"  (D.T.  Landin,  C.W.  Macosko)  submitted  to 
Macromolecules,  1986. 

'Structural  and  Rheological  Changes  During  Epoxy  Crosslinking"  (S.A.  Bidstrup  and  C.W 
Macosko),  to  appear  in  proceedings  from  symposium  on  Crosslinked  Epoxies,  Prague  July  14-17, 
1986. 

Previously  Reported 


"Molecular  Weight  Relations  for  Crosslinking  of  Chains  with  Length  and  Site  Distribution"  (C.W 
Macosko  and  D.R.  Miller),  Journal  of  Polymer  Science,  Physics  Edition,  accepted  for  publication 
1986. 

"Network  Parameters  for  Crosslinking  of  Chains  with  Length  and  Site  Distribution."  (C.W. 
Macosko  and  D.R.  Miller),  Journal  of  Polymer  Science,  Physics  Edition,  accepted  for  publication 
1986. 

"Structure-Rheology  Relations  for  Model  Epoxy  Networks"  (Sue  Ann  Bidstrup  and  C.W. 
Macosko)  presented  at  SAMPE  meeting,  April  7-10,  1986. 

"Rheological  Changes  during  Crosslinking"  (C.W.  Macosko),  British  Polymer  Journal,  JLZ, 
239-245  (1985). 

"Heat  Transfer  and  Property  Development  in  Liquid  Silicone  Rubber  Molding"  (C.W.  Macosko 
and  L.J.  Lee)  Rubber  Chem.  Tech.,  1£,  436-448  (1985). 

"Viscosity  Rise  During  Free  Radical  Crosslinking  Polymerization  With  Inhibition"  (V.M.  Gonzalez 
and  C.W.  Macosko)  Journal  of  Rheology,  2JL  259-272  (1985). 

"Ring-chain  Competition  Kinetic  Models  for  Linear  and  Nonlinear  Step-reaction  Copolymerizations 
"(C.  Sarmoria,  E.M.  Vallds  and  D.R.  Miller)  Makromol.  Symp.,  2,  69-87  (1986). 

"Rheological  Changed  During  the  Copolymerization  of  Vinyl  and  Divinyl  Monomers"  (D.T. 

Landin  and  C.W.  Macosko)  in  "Characterization  of  Highly  Crosslinked  Polymers,"  eds.  S.S. 
Labana  andR.A.  Dickie,  ACS  Symposium  Series,  243.  33-46  (1984)  Washington,  D.C. 

"Small  Strain  Modulus  of  Model  Trifunctional  Polydimethyl  Siloxane  Networks"  (E.M.  Valles, 

E.J.  Rost  and  C.W.  Macosko)  Rubber  Chem.  Tech.  12,  55-62  (1984). 

"Structural  and  Rheological  Changes  During  Epoxy  Crosslinking"  (S.A.  Bidstrup  and  C.W. 
Macosko)  Soc.  Plast.  Eng.  Tech.  Papers,  IQ,  278-281  (1984). 


"Molecular  Weight  Change  During  Vinyl-Divinyl  Copolymerization  (D.T.  Landin  and  C.W. 
Macosko;  presented  at  European  Networks  Group  Meeting,  Manchester,  England,  September. 
1984. 
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Calculation  of  Average  Network  Parameters  Using  Combined  Kinetic  and  Markovian  Analysis" 
(D.R.  Miller  and  C.W.  Macosko),  submitted  for  publication  in  "Biological  and  Synthetic 
Networks",  Elsevier,  to  appear  1987. 

Cvclization  and  Reduced  Reactivity  of  Pendant  Vinyls  During  the  Copolymerization  of  Methyl 
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"Structural  and  Rheological  Changes  During  Epoxy  Crosslinking"  (S.A.  Bidstrup  and  C.W. 
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Calculation  of  Average  Network  Parameters 
Using  Combined  Kinetic  and  Markovian  Analysis 


D.  R.  Miller,  Department  of  Operations  Research,  George  Washington  University, 
Washington  D.C.  20052 

C.  W.  Macosko,  Department  of  Chemical  Engineering  and  Materials  Science,  University  of 
Minnesota,  Minneapolis,  MN  55455 


Theoretical  models  are  useful  for  calculating  various  average  parameters  during 
network  formation  such  as  weight  average  molecular  weight  before  gelation,  the  gel  point, 
and  weight  fraction  solubles,  and  concentration  of  pendant  chains  and  elastically  effective 
junctions  after  gelation.  These  models  and  their  methods  of  solution  may  take  various  forms; 
the  most  common  are  kinetic  models  and  Markovian  models.  The  kinetic  approach  applies  to 
a  larger  class  of  systems;  however,  it  may  become  intractable,  especially  for  branching 
polymerizations,  because  a  large  number  of  differential  equations  must  be  solved.  The 
Markovian  approach  (sometimes  referred  to  as  the  statistical  or  cascade  approach)  is  attractive 
because  it  permits  more  direct  calculation  of  average  properties  by  exploiting  the  recursive 
structure  of  the  material.  An  optimal  approach  would  combine  the  two  methods.  This  paper 
surveys  some  systems  which  can  be  efficiently  analyzed  with  a  combined  methodology. 

The  kinetic  approach  and  the  Markovian  approach  can  be  contrasted  by  considering 
the  gel  point  of  an  ideal,  A^  stepwise  homopolymerization.  Before  gelation,  the  kinetic 
approach  models  the  concentration  of  each  size  of  molecule.  If  Pi  represents  a  polymer 
consisting  of  i  Af  monomers  then: 


P  +  P 

‘i  j 


at  a  rate  proportional,  ideally,  to  [(fi-i-t-l)(fj-j+l)]  for  i,  j  >  1.  Differential  equations  are 
solved  to  find  the  concentration  of  different  molecular  sizes  as  a  function  of  time  or 


conversion.  Analytical  solution  of  these  equations  is  quite  difficult  for  even  the  simplest 
systems  (Dusek,  1985a;  Ziff,  1980;  Ziff  and  Stell,  1980)  and  appears  to  be  impossible  for 
complex  systems,  and  thus  the  equations  must  be  solved  numerically.  From  the 
concentration  of  different  species  the  weight  average  molecular  weight  is  computed. 


Gelation  occurs  when  the  weight  average  molecular  weight  diverges. 

The  Markovian  approach  is  based  on  the  fact  that  the  monomers  combine 
randomly  in  an  ideal  homopolymerization,  i.e.  all  intermolecular  AA  bonds  are  equally  likely 
to  form.  This  results  in  a  branching  structure  which  has  a  recursive  structure  with 
lack-of-memory  which  can  be  exploited  by  various  analysis  methods  (Flory,  1953;  Gordon, 
1962;  Macosko  and  Miller,  1976;  Burchard,  1983).  For  illustration  we  use  the  "in-out  ' 
recursion:  let  EWAout  equal  the  expected  weight  seen  looking  out  from  a  randomly  chosen 
A-site  (directions  1  in  Figure  1).  At  conversion  a: 

EWAout  =  oEWAin  +  (1-  a)  0  (2) 

where  EWAm  is  the  expected  weight  looking  into  a  random  A-site  (directions  2  in  Figure  1). 

EwAin  =  Ma+  (f-1)  EWAout  (3) 

Finally  the  weight  average  molecular  weight  is: 

Mw  =  MA  +  fEWA0Ut  (4> 


Solving  the  above  three  equations  gives: 


MW“MA  +  aMA/l  Wf*  D]  (5) 

and  gelation  occurs  at  ot^.  =  1  /  (f  -1),  which  is  identical  to  the  answer  computed  using  kinetic 
analysis. 

The  Markovian  analysis  assumes  that  the  species  combine  in  a  totally  random  way. 
There  are  many  systems  for  which  this  is  not  true.  In  these  cases  we  recommend  combining 
the  two  methods.  This  idea  is  not  new;  several  authors  have  used  it  implicitly  and  Dusek 
( 1986)  has  explicitly  presented  it  as  a  formal  approach  to  analyzing  different  systems. 

In  a  combined  analysis  the  aspects  of  the  reaction  which  are  not  totally  random  are 
treated  kinetically  to  obtain  "superspecies"  which  are  combined  randomly  and  thus  can  be 
treated  using  a  Markovian  analysis; 

Initial  Species  — ►  Superspecies  -*  Product 

nonrandom  random 

combinauon  combination 


The  two  steps  in  this  model  may  either  correspond  to  two  distinct  sequential  steps  or  they 


may  be  occurring  simultaneously,  in  which  case  the  construct  is  merely  a  convenient  way  to 
think  about  the  process.  This  paper  discusses  several  different  systems  whose  average 
properties  can  be  computed  efficiently  using  a  combined  kinetic -Markovian  approach. 


CoDQlvmerization  with  Substitution  Effect 

Consider  an  A3  plus  B2  copolymerization  for  which  there  is  a  first  shell  substitution 
effect  on  the  A3  monomers.  This  system  has  been  analyzed  by  numerous  authors,  including 
Miller  and  Macosko  (1980).  By  explicitly  recognizing  the  combined  kinetic-Markovian 
approach  some  clarity  is  revealed;  it  is  a  good  example  for  illustrating  the  concepts. 

We  define  seven  superspecies:  four  types  of  A-monomer  and  three  types  of 
B- monomer.  Let: 

A3‘  *  A-monomer  with  i  reacted  sites,  i  **  0, 1, 2, 3 

B2‘  *  B-  monomer  with  i  reacted  sites,  i  =  0, 1, 2 

kj  =  rate  constant  for  reaction  of  A3‘  functionality  with  a  B  functionality,  i  =  0,1,2. 
From  a  kinetic  analysis  we  can  find  the  concentration  of  each  of  the  seven  superspecies  by 


solving  the  following  system  of  kinetic  equations: 

1  A3°  (t)  =  -3kg  A30  (t)  (2B2°(t)  +  B21  (t))  (6) 

1 A3!  (t)  =  (3kg  A3°  (t)  -  2k1A31  (t))(2B2°(t)  +  B2l(t))  (7) 

A32  (t)  =  (  2k1A31  (t)  -  k^2  (t))(2B2°(t)  4-  B2*  (t))  (8) 

~  A33  (t)  =  k2A32  (t)  (2B2°(t)  +  B2l  (t))  (9) 

-|t  B2°  (t)  *  -  2B2°(t)  (3kg  A3°  (t)  +  2k1A3I  (t)  +  k2A32  (t))  ( 10) 

J-  B21  (t) » ( 2B2°(t)  -  B21  (t))  (3kg  A3°  (t)  +  2k!  A31  (t)  +  k^2  (t))  (11) 

J-  B22  (t)  =  B2l  (t)  (3kg  A3°  (t)  +  2klA31  (t)  +  l^Aj2  (t))  (12) 


The  seven  superspecies  are  shown  in  Figure  2.  The  combination  is  random:  any  pair 
of  reacted  A-site  and  reacted  B-site  is  equally  likely  to  combine.  Thus  instead  of  thinking  of 
the  system  as  two  species,  A3  and  B-,,  combining  together  in  a  rather  complicated,  not 


completely  random  way  (substitution  effect),  we  think  of  seven  superspecies  which  combine 
in  a  completely  random  way.  Furthermore  the  seven  superspecies  go  to  complete  reaction. 
The  resulting  structure  is  a  branching  process  with  seven  types  of  components:  it  can  be 
analyzed  using  any  of  several  methods  (the  pgf  approach  of  Gordon  (1962),  the  recursive 
approach  of  Macosko  and  Miller  (1976),  the  analysis  of  generation  sizes  as  mentioned  by 
Burchard  (1983),  or  even  the  Flory-Stockmayer  combinatorial  approach  (Rory,  1953)).  The 
point  is  that  by  the  construct  of  superspecies  it  is  possible  to  identify  a  Markovian  branching 
process  which  can  then  be  analyzed  using  existing  techniques. 


Homopolvmerization  with  Substitution  Effect 


Consider  an  Af  homopolymerization  with  first-shell  substitution  effect  At  first 
glance  this  might  be  considered  as  a  special  case  of  the  above  copolymerization  example,  but 
in  general  it  is  more  complicated.  For  this  case  it  is  necessary  to  define  more  complicated, 
superspecies  of  A- monomers,  if  our  goal  is  a  Markovian  branching  structure. 

For  a  general  first  shell  substitution  effect  homopolymerization  we  have  a  matrix  of 
kinetic  rates: 


^"*00 


02 


K  = 


10 


Ml 


M2 


^.f-l 

k. 


M,  f-1 


(13) 


L  kf-l.O  kf-l,  1  kf-l,  2‘  •  -kf-l.f-l  -I 


where 


kjt  j  »  rate  at  which  a  monomer  with  i  reacted  sites  forms  a  bond  with  a  monomer 
with  j  reacted  sites 


and 


Jr  .  =  k  •  • 

i.J  J>i* 


We  define  the  superspecies  by  labelling  each  reacted  site  of  a  monomer  with  a  pair  of  integers 
(m,n),  where  m  signifies  that  it  is  the  m*  site  to  react  on  the  monomer  and  n  signifies  that  the 
site  with  which  it  reacts  is  the  n^  site  to  react  on  the  monomer  to  which  it  belongs,  l£m,n£f. 
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The  superspecies  for  an  A2  homopolymerization  with  first  shell  substitution  effect  are  shown 
in  figure  3;  there  are  seven  of  them.  The  situation  is  similar  for  general  f  except  the  number 
of  superspecies  increases  rapidly:  for  f  =  3  there  are  40  superspecies. 

These  superspecies  combine  randomly  into  a  Markovian  branching  structure:  an  (m, 
n)  site  combines  with  an  (n,  m)  site.  By  having  such  complicated  superspecies  we  have 
created  a  system  without  longer  range  dependencies  which  exist  in  other  models  of 
homopolymerization  with  general  first-shell  effect 


Consider  a  system  which  forms  by  crosslinking  of  premade  chains  which  have 
arbitrary  distributions  of  mass  and  functionality.  In  this  case  the  chains  play  the  role  of 
supeispecies  so  the  two  stages  are  distinct:  manufacture  of  the  functional  polymer  (typically 
kinetically  controlled)  and  then  random  crosslinking  using  radiation  or  a  small  monomer 
(typically  Markovian).  This  case  is  treated  by  Miller  and  Macosko  (1986a,  1986b). 


Consider  an  A3  +  B2  copolymerization  in  which  intramolecular  bonds  are  allowed  to 
form.  Sarmoria  et  al.  (1986a)  have  analyzed  a  special  case  of  this  system  in  which  the  only 
allowable  intramolecular  bonds  are  those  which  close  a  cycle  consisting  of  a  single  B2 
monomer  forming  two  bonds  with  a  single  A^  monomer.  They  define  13  superspecies 
which  correspond  to  the  possible  configurations  within  one  shell  of  A-monomers  and  four 
superspecies  corresponding  to  B -monomers.  The  concentration  of  the  various  superspecies 
are  kinetically  modelled  assuming  a  Gaussian  end-to-end  distance  and  competing  ring-chain 
reactions.  The  superspecies  are  then  combined  randomly  into  a  Markovian  branching 
structure. 

Sarmoria  et  al.  (1986b)  have  extended  the  above  model  to  include  larger  cycles. 
Cycles  involving  a  single  B2  are  modelled  as  before,  but  larger  cycles  are  modelled  using  the 
spanning  tree  approximation  of  Gordon  and  Scantlebury  (1966).  In  this  case,  as  the 
polymerization  is  proceeding,  intermolecular  bonds  are  labelled  with  a's,  intramolecular 
bonds  are  labelled  with  o's,  and  the  unreacted  sites  with  co’s.  The  intramolecular  bonds  on 
cycles  with  more  than  one  B2  are  then  broken  (each  site  retaining  its  o-label).  This  results  in 
a  tree  structure  made  up  of  22  different  superspecies.  As  before  the  superspecies  reflect  all 
possible  configurations  within  the  first  shell  surrounding  an  A-monomer.  The  22 


superspecies  are  shown  in  Figure  4.  Numerous  variations  on  these  models  are  possible. 


Step  Reactions  which  Produce  a  Secondary  Site 

Consider  an  Af  +  Bg  copolymerization  such  as  an  f-functional  epoxy  and  a 
g-functional  cresol;  such  a  system  is  used  for  transitor  encapsulation.  In  these  systems  bond 
formation  between  two  sites  creates  a  new  site.  One  of  the  initial  reactants  can  react  with  this 
site,  often  in  a  chainwise  manner.  The  block  of  connected  sites  generated  from  one  initial 
reaction  is  one  of  the  superspecies,  although  in  this  context  it  might  be  more  appropriate  to 
call  them  "super-bonds."  Figure  5  shows  the  situation. 

Bokare  and  Gandhi  (1980)  and  Dusek  (1985b)  have  analyzed  the  reaction  of  epoxy 
and  diamine  with  polyetherification  using  a  combined  kinetic-Markovian  approach.  Riccardi 
and  Williams  (1986  a,b)  have  also  studied  these  systems.  Other  such  systems  include 
isocyanate  and  hydroxyl  with  allophanate  formation. 


Anionic  Copolvmerization  of  Vinvl  with  Divinvl  Monomers 

Worsfold  (1970)  and  Lutz  et  al.  (1982)  describe  anionic  copolymerizations.  Lutz  et 
al.  consider  a-methylstyrene  (a-MS)  plus  1,3-diisopropenylbenzene  (1,3-DIB).  This 
system  is  depicted  on  Figure  6;  the  chemical  structure  is  shown  only  in  unreacted  species. 
The  superspecies  are  the  living  polymer  chains,  when  a  1,3-DIB  monomer  is  part  of  a  chain 
its  other  group  may  be  unreacted  (labelled  with  co  in  Figure  6)  or  reacted  in  which  case  the 
monomer  acts  as  a  crosslink  (and  the  other  group  is  labelled  with  an  a).  The  relevant 
properties  of  the  superspecies  are  computed  from  a  kinetic  analysis;  these  properties  arc  very 
similar  to  those  needed  in  the  analysis  of  the  crosslinking  of  ready-made  chains  studied  by 
Macosko  and  Miller  (1986a,  1986b).  The  superspecies  are  then  randomly  combined  together 
into  a  Markovian  structure  and  the  product  analyzed  by  branching  theory.  Dusek  and 
Somvarsky  (1985)  have  also  studied  similar  systems  using  a  combined  kinetic-Markovian 
analysis. 
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figure  captions 


Figure  1.  Schematic  of  a  typical  branched  molecule  after  degree  of  polymerization  a. 


Figure  2.  Initial  species,  superspecies  and  product  for  A3  +  B2  copolymerization  with  first 
shell  substitution  effect  on  A3. 


Figure  3.  Species,  superspecies  and  product  for  A2  homopolymerization  with  first  shell 
substitution  effect 


Figure  4.  Superspecies  for  A3  +  B2  copolymerization  with  cycles. 


Figure  5.  Af  +  Bg  copolymerization  with  secondary  B-site  produced. 


Figure  6.  a-MS  +  1, 3-DEB  copolymerization. 
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Cyclization  and  Reduced  Reactivity  of  Pendant  Vinyls  During 
the  Copoi ymeri zati on  of  Methyl  Methacrylate  and  Ethylene 

Glycol  Dimethacrylate 

D.  T.  Landin  and  C.  W.  Macosko 


Abstract 

Methyl  methacrylate  was  copolymerized  with  small 
amounts  of  ethylene  glycol  diraethacrylate .  Monomer  and 
pendant  vinyl  conversion  as  a  function  of  time  were  measured 
up  to  the  gel  point.  Pendant  vinyl  conversion  was 
.etermined  by  1 H  NMR .  Plots  of  pendant  conversion  vs. 
monomer  conversion  exhibit  a  positive  v-intercept  indicating 
the  tendency  to  cycle  during  the  formation  of  a  primary 
chain.  For  bulk  systems,  this  amounted  to  approximately 
3-  4%  o'  pendant  vinyls  and  the  proportion  in  cycles 
increased  wi th  dilution.  The  slope  of  the  pendant  vs . 
monomer  conversion  plot  is  attributed  to  the  formation  of 
crosslinks  and  subsequent  cycles.  A  kinetic  model  is 


INTRODUCTION 


Understanding  the  behavior  of  pendant  vinyls  is  key  to 
the  study  of  network  formation  by  addition  polymerization 
involving  multifunctional  vinyl  monomers.  During  such  a 
polymerization,  pendant  vinyls  are  created  when  the  first 
vinyl  unit  of  a  multifunctional  vinyl  monomer  reacts  and  the 
group  is  added  to  a  polymer  chain.  The  pendant  vinyl  can 
then  react  or  remain  pendant.  The  reactivity  of  the  pendant 
vinyl  may  be  the  same  as  that  of  the  monomeric  vinyl  or  it 
may  be  increased  or  reduced.  If  it  reacts  with  a  growing 
chain  to  which  it  is  not  already  chemically  attached,  it 
forms  a  crosslink.  Cyclization  occurs  if  the  pendant  vinyl 
is  already  attached  to  the  chain  with  which  it  reacts. 

This  behavior  of  pendant  vinyls  affects  the  prediction 
of  network  properties  such  as  crosslink  density,  sol 
fraction,  and  critical  conversion  for  gelation.  The  first 
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theories  which  connected  structural  properties  and  extent  of 
reaction  were  those  of  Flory1'*  and  Stockmayer 3  *  4 .  These 
theories  assumed  equal  reactivity  of  monomeric  and  pendant 
vinyls  and  all  reacted  pendant  vinyls  were  crosslinks  (no 
cvc 1 i za t i on ) .  Walling5  showed  that  for  the  methyl 
methacrylate  -  ethylene  glycol  dimethacrylate  and  vinyl 
acetate  -  divinyl  adipate  systems  there  was  a  large 
discrepancy  between  the  ge]  point  calculated  from  the 
FI ory-Stockmayer  theory  and  that  found  experimentally.  He 
attributed  the  discrepancy  to  diffusion  control  of  the 


reaction. 


In  a  sci  ics  of  papers,  Gordon  and  Roe 6  disputed 


this  explanation  and  instead  attributed  the  discrepancy  to 
the  formation  of  cycles.  There  is  extensive  experimental 
evidence  of  the  inadequacy  of  the  FI ory-Stockmayer  theory 
for  several  systems'-15.  Two  explanations  for  the 
discrepancy  between  observed  and  theoretical  gel  points  are 
the  formation  of  cycles16  and  a  reduced  reactivity  of  the 
pendant  functional  group.14*  17 

In  this  work,  conversion  of  monomeric  and  pendant 
vinyls  and  the  gelpoint  are  determined  experimentally  for 
the  vinyl-di vinyl  system  of  methyl  methacrylate  with  small 
amounts  of  ethylene  glycol  dimethacrylate.  A  simple  model 
is  proposed  which  includes  the  possibility  for  both 
cyclization  and  reduced  reactivity.  This  provides  a 
framework  within  which  to  organize  and  study  the  data. 

EXPERIMENTAL 

4 

Methyl  methacrylate  ( MMA )  and  ethylene  glycol 
dimethacrylate  ( EGDMA )  were  obtained  from  Aldrich.  MMA  was 
washed  with  a  10%  aqueous  potassium  hydroxide  solution  and 
EGDMA  was  used  as  received.  The  monomers  were  mixed  to  give 
the  chemical  systems  shown  in  Table  1. 

Azob i s i sobu t y ron i t r i 1 e  (A1BN)  from  Kodak  was  added  to  the 
solutions  in  the  amount  0.3%  by  weight.  To  some  solutions, 
1%  by  volume  dodecaneth i o 1  was  added  as  a  transfer  agent. 

The  reaction  mixture  was  poured  into  several  6  inch  glass 
test  tubes  each  capped  with  a  rubber  septum.  The  tubes  were 
placed  in  a  60°C  water  bath  for  varying  amounts  of  time. 


Tubes  were  removed  from  the  bath  and  polymer  was 
precipitated  in  methanol  to  which  a  small  amount  of  dipheny 
pjcryl  hydrazyl  (DPPH)  had  been  added. 

The  conversion  of  monomer  was  determined  by  drying  and 
weighing  the  precipitated  polymer.  Conversion  vs.  time  data 
for  the  polymerization  under  these  conditions  are  shown  in 
Figure  1.  Neither  the  addition  of  EGDMA  or  the  addition  of 
transfer  agent  changed  the  conversion  vs.  time  results.14 

Gel  points  were  determined  by  observing  the  time  at 
which  the  reaction  mixture  would  no  longer  flow.  The  gel 
point  conversion  was  then  obtained  from  the  conversion  vs. 
time  plot.  These  results  are  listed  in  Table  2. 

Proton  nuclear  magnetic  resonance,  >H  NMR ,  was  used  to 
determine  the  conversion  of  pendant  vinyls.  The 
precipitated  polymer  was  dissolved  in  chloroform  and  then 
reprecipitated  in  methanol.  The  polymer  was  dried  and  then 
dissolved  in  deuterated  chloroform  with  1%  v/v 
tetramethyl s i lane  (TMS).  This  solution  was  used  for  NMR 
analysis.  NMR  spectra  were  obtained  with  a  Nicolet  NT-300 
WB  (Nicolet  Analytical  Instruments).  A  one-pulse  sequence 
was  used  with  a  45°  pulse,  2  second  delay  time,  and  2.05 
second  acquisition  time.  Approximately  100  acquisitions 
were  obtained  per  spectrum.  Quadrature  detection  was  used 
and  the  spectral  window  was  4000  Hz.  The  spectra  were 
machine  integrated. 

An  NMR  spectrum  for  sample  2N  is  shown  in  Figure  2.  A 
blowup  of  the  region  from  4.0  to  6.5  ppm  is  shown  in  Figure 


3.  The  3  peaks  from  0.7  to  1.25  ppm  are  from  methyl  groups 
adjacent  to  reacted  vinyl  groups.  The  peaks  from  -5.0  to  4.5 
ppm  are  from  the  -OCH 2-  protons  on  EGDMA .  The  two  peaks  at 
5.65  and  6.2  ppm  are  from  unreacted  vinyl  groups.  The  peak 
at  3.55-3.7  ppm  is  due  to  -OCH 3  protons  from  MMA . 

The  areas  under  the  -OCH2-  and  -OCH 3  peaks  are  used  to 
determine  copolymer  composition  through  the  equation 


leobMA]  ^  I!  a  / 

[MMA]  At-OUj^/S  ‘ 

where  A(-OCHa-)  and  A(-0CH2)  are  the  areas  under  the  peaks. 
Using  Equation  1,  the  polymer  ratio  is  found  to  be  identical 
to  the  monomer  ratio  (Table  1)  indicating  monomeric  vinyl 
reactivities  are  equal  as  found  by  previous  workers15. 
Pendant  conversion  is  given  by  the  equation 


>,  _  X  Af-cgcujt/g 

Ya  (-OC Hy)7T" 

where  pp  denotes  pendant  conversion. 

t 

In  some  cases,  it  was  very  difficult  to  eliminate  all 
monomer  from  polymer.  When  this  happened,  the  spectrum 
appeared  as  in  Figure  4.  It  was  difficult  to  separate  the 
monomer  and  polymer  -OCH2-  peaks  so  an  alternative  analysis 
was  used.  This  analysis  used  the  equation 


*  ‘  —ZTfi -  £) 
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/S  *  2  A(-  C:cw. )/ 2- 


Intrinsic  viscosity  was  used  to  determine  the  molecular 
weight  of  MMA  with  no  HJDMA  polymerized  under  the  same 


f'nnd  i  t  i  ons  as  those  with  f  GDMA  . 


Meas'ii  erne n t  h  were  made  at 


3(1 0 C  with  a  0.46  mm  diameter  Ubbelhode  viscometer.  Ethyl 
acetate  (EtO Ac)  was  used  as  solvent.  .  Using  the  equation 

t7!]  =  KM* 

where  K  and  a  for  MMA  in  EtOAc  at  30 °C  are1  ‘ 

K  r  x  ic> 

ol  = 

The  molecular  weights  were 

M  =  2.47x10s 

for  MMA  without  any  transfer  agent  and 

M  =  5.21x10* 

for  MMA  with  transfer  agent. 

RESULTS  AND  DISCUSSION 

Plots  of  pendant  vinyl  conversion  vs.  monomer 
conversion  are  shown  in  Figures  5-7.  Figure  5  shows  the 
effect  of  increasing  amounts  of  EGDMA.  Lines  have  been 
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drawn  through  the  data  points  using  a  least  squares  fit  to 
better  show  the  trend  of  the  data.  All  systems  show  a 
positive  y-intercept.  This  indicates  that  in  the  limit  of 
zero  conversion,  pendant  vinyls  have  reacted  or  in  other 
words  as  the  first  polymer  chains  are  forming,  pendant 
vinyls  are  reacting.  Since  crosslinking  is  a  second  order 
reaction  and  in  the  limit,  of  zero  conversion  polymer  chains 
are  infinitely  dilute,  this  loss  of  pendant  vinyls  must  be 


due  to  cyclwat ion.  Cycl ization  during  the  formation  of  a 
primary  chain  will  be  called  primary  cyclization . 

From  Figure  5,  it  can  be  seen  that  within  this  range  of 
crosslinker  concentration ,  increasing  the  EGDMA  does  not 
change  the  proportion  used  up  in  primary  cycles.  The  effect 
of  a  different  primary  chain  length  is  shown  in  Figure  6. 
Samples  with  and  without  transfer  agent  have  the  same 
proportion  of  pendant  vinyls  in  primary  cycles,  however  the 
slopes  are  much  different.  Figure  7  shows  the  effect  of 
dilution.  A  greater  proportion  of  pendant  vinyls  are  used 
in  primary  cycles  as  the  dilution  is  increased. 

The  y-intercept  and  slope  from  each  pendant  conversion 
curve  and  the  gel  point  conversion  for  each  chemical  system 
investigated  are  listed  in  Table  2. 

In  developing  a  model  to  describe  this  behavior,  the 
slope  of  the  pendant  conversion  curve  is  attributed  to  the 
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formation  of  crosslinks  and  to  secondary  cycles.  Secondary 
cycles  are  those  formed  with  a  pendant  vinyl  attached  to  a 
polymer  chain  which  has  just  been  connected  to  the 
propagating  radical  chain  by  a  crosslink.  The  process  of 
crosslinking  and  secondary  cyclization  is  illustrated  in 
Figure  8.  The  contributions  of  crosslinking  and  secondary 
cyclization  to  the  slope  can  be  separated  at  the  gel  point 
since  the  conversion  of  pendant  vinyls  to  crosslinks 
necessary  for  gelation  is  known  theoretically. 
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A  detailed  d  e  \  el  i  >  pnu*  ri  1  of  the  kinetic  mode]  is  found  in  the 
Appendix.  Three  constants  oi  concern  are: 


r  i  -  react i  v  i  t  >  ratio  of  pendant  vinyl  to  monomeric  vinyl 

kPP  -  constant  for  primary  cyclization 

kpS  -  constant  for  secondary  cyclization. 

The  primary  cyclization  constant  is  obtained  from  the 
y-intercept  of  the  pendant  conversion  vs.  monomeric 
conversion  plot.  The  slope  of  this  plot  is  not  enough 
information  to  determine  n  and  kpi  alone.  A  further  piece 
of  information  is  needed.  This  is  obtained  by  using  network 
theory.  The  pendant  conversion  to  crosslinks  necessary  for 
gelation  is  obtained  theoretically  from  the  equation20 

l 

where  p  is  the  conversion  of  pendant  vinyls  to  crosslinks 
necessary  for  gelation,  q  is  the  probability  that  a 
monomeric  unit  incorporated  in  a  polymer  chain  is  not  a 
chain  end  and  a«  is  the  proportion  of  vinyl  functional 
groups  which  are  on  divinyl  monomer  units. 

To  obtain  a  unique  pair  {kp«,rj)  of  values  which  gives 
both  the  correct  slope  and  the  correct  crosslink  conversion 
at  gelation,  kPI  s.  r  j  giving  the  correct  slope  and  kP.  vs. 
ri  giving  the  crosslink  conversion  are  plotted. 
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The  kps  vs.  r3  plot  for  pairs  giving  the  correct,  slope 
of  the  pendant  conversion  plot  is  determined  by  using  the 
kinetic  model  developed  in  the  Appendix.  A  value  for  r3  is 
chosen  and  the  slopes  obtained  using  three  different  values 
for  kps  with  the  r3  value  are  used  to  find  the  kps  value 
giving  the  correct  slope  by  linear  regression.  This  is 
repeated  for  several  different  r3  values. 

The  kps  vs.  r3  plot  for  pairs  giving  the  correct 
conversion  to  crosslinks  at  gelation  is  also  determined 
using  the  kinetic  model.  Once  again,  a  value  for  r3  is 
chosen.  The  conversion  to  crosslinks  at  the  experimentally 
determined  monomer  conversion  of  gelation  is  calculated  with 
this  r3  value  and  three  different  kp(  values.  The  kp(  value 
giving  the  crosslink  conversion  predicted  by  Equation  4  is 
determined  by  linear  regression.  This  is  repeated  for 
several  r3  values. 

The  intercept  of  the  two  curves  gives  the  values  for  r3 
and  kPi.  Figure  9  shows  this  plot  for  sample  2T.  With  and 
without  transfer  agent,  the  reactivities  are  similar  but  the 
amounts  of  secondary  cyclization  are  very  different.  Table 
2  shows  r3  and  kp«  as  determined  for  several  of  the  systems 
investigated.  When  the  amount  of  crosslinker  is  held 
constant  and  the  chain  length  varied,  the  pendant  reactivity 
is  fairly  constant  while  the  number  of  secondary  cycles  per 
crosslink  increases  with  increasing  chain  length.  This  may 
be  so  because  a  longer  chain  has  a  greater  concentration  of 
pendant  vinyls  surrounding  a  radical  end.  In  the  case  of 


sys  terns 


2 S'  Find  2T,  the  system  with  transfer  agent  has  a  kP« 
equal  to  2.7  and  the  system  without  transfer  agent  a  value 
of  19.7.  This  is  reasonable  because  the  primary  chain 
length  in  the  system  without  transfer  agent  is  nearly  six- 
times  larger  than  in  the  system  with  transfer  ager  Both 
the  pendant  reactivity  and  the  number  of  secondary  cycles 
change  as  the  primary  chain  length  is  held  constant  and  the 
amount  of  crosslinker  is  varied.  Pendant  reactivity 
decreases  as  the  amount  of  cosslinker  is  increased.  This 
is  probably  due  to  an  increase  in  steric  screening  effects 
caused  by  increasing  the  number  of  branches  in  the  vicinity 
of  a  pendant  vinyl  as  well  as  decreasing  the  mobility  of  the 
macromolecule.  The  number  of  secondary  cycles  increases  as 
the  amount  of  crosslinker  is  increased.  This  is  once  again 
due  to  the  greater  concentration  of  pendant  vinyls 
surrounding  a  radical  end. 

Once  the  constants  have  been  determined,  the  model  can 
be  used  to  predict  total  conversion  of  pendant  vinyls 
and£ onve rs ion  to  crosslinks,  primary  cycles,  and  secondary 
cycles.  These  predicted  conversions  are  shown  as  a  function 
of  monomer  conversion  in  Figure  10  for  system  2T.  These 
conversions  can  now  be  used  in  a  branching  theory  to  predict 
molecular  parameters. 1  s  >  21 

CONCLUSIONS 

Vinyl -di vi nyl  network  formation  has  been  investigated 
bv  determining  the  conversion  of  pendant  vinyls  as  a 


function  of  monomeric  \  i nyl  conversion,  using  this 
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information  to  shape  the  concept  of  cyclization  and 
crosslinking,  and  formal  izing  these  concepts  by 
incorporating  them  in  a  kinetic  model  for  the  rate  of 
polymerization.  Plots  of  pendant  conversion  vs.  monomer 
conversion  exhibited  n  positive  y-intercept  indicating  the 
loss  of  pendant  vinyls  during  the  formation  of  a  primary 
chain.  This  loss  was  attributed  to  the  formation  of  primary 
cycles.  Primary  cyclization  increased  with  increasing 
dilution.  The  slope  of  the  plot  indicated  the  formation  of 
secondary  cycles  with  the  f orma tion  of  crosslinks.  The 
kinetic  model  along  with  experimental  conversion  curve  data 
provided  a  means  to  estimate  the  number  of  secondary  cycles 
and  the  reactivity  of  the  pendant  vinyl.  The  amount  of 
secondary  cyclization  increased  with  increasing  length  of 
the  primary  chain.  The  pendant  vinyl  reactivity  was  found 
to  be  approximately  half  that  of  monomeric  vinyls. 

This  approach  to  understanding  network  formation  is 
that  of  developing  a  model  which  is  consistent  with  the 
available  experimental  data,  and  as  is  true  with  models,  it 
is  only  an  approximation.  Primary  cyclization  has  been 
assumed  to  occur  at  a  constant  rate.  This  is  probably  not 
so  since  as  the  reaction  progresses,  the  concentration  of 
pendant  vinyls  not  connected  by  a  primary  chain  to  a 
particular  radical  end  increases  arid  provides  competition. 
Pendant  reactivity  in  this  model  is  not  that  of  an  isolated 
vinyl  in  a  very  dilute  system.  Instead  it  must  include 
contributions  from  many  different  sources  in  one  term.  It 


t . 
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lb  not  surprising  then  that  it  is  less  than  the  reactivity 


of  monomeric  vinyls. 
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APPEND]  X 


With  the  concept  of  primary  and  secondary  cyclization 
in  mind,  a  kinetic  scheme  for  the  vi ny 1 -d i v iny 1 
copolymerization  can  be  written.  In  this  scheme, 

A  =  vinyl  on  mono vinyl  monomer 
B  =  vinyl  on  divinyl  monomer 
P  =  pendant  vinyl 

a  =  vinyl  in  polymer  from  monovinyl  monomer 
b  =  first  vinyl  from  divinyl  monomer  incorporated  in 
pol ymer 

p  =  pendant  vinyl  incorporated  in  polymer 
pp  =  pendant  vinyl  in  primary  cycle 
px  =  pendant  vinyl  in  crosslink 
ps  =  pendant  vinyl  in  secondary  cycle 


The  propagation  steps  are  written  as 

A*  4  A  — »  A*  (kn) 

A*  4  "B  — S*  (  ki^ 

*  P  -►  P* 

B*  ♦  A  A*  (  k.1,) 


The  equHt  if.ii.1-  fi'f  the  tame  rate  of  appearance  and 
d a sappea ranee  of  species  A  and  B  are 
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-  *  a fk  1-4*1161  tjulP-lW). 

Pendant  vinyls  are  both  produced  and  consumed  so  an  equation 
for  the  time  rate  of  change  has  terms  for  both. 

Contributions  to  the  term  for  disappearance  come  from 
reaction  to  crosslinks,  primary  cycles,  and  secondary 
cycles.  The  equation  for  P  becomes 

■  *  0  *  *  Mb-Hp]  ■> h5[p']lp3 

-  ¥  5T 

where  kPP  is  the  fraction  of  pendant  vinyls  in  primary 
cycles  and  kP,  is  the  average  number  of  secondary  cycles 
formed  per  crosslink.  Differential  equations  are  also 

J 

written  for  the  appearance  of  species  in  the  polymer.  These 
are 

_  ^3  c*) 

d  t  dt- 
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di 
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Combining  Equations  5  and  6  gives 
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Assuming  quasi-steady  state  and  using  the  accompanying 
equat i ons 

kjA-UE]  -  k1(  IS*1U] 

kn  Ip*]  [a]  *  k,sU*][pj 


Equation  11  becomes 

5^^  _  U1  ( kii  /  ki;  La]  +  s  +  ^'i/ ku  \p~l ) 
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If  r,3  is  defined  as 


Equation  15  becomes 
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Equations  G  and  7  combine  to  give 
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Expressions  lor  the  relatrn  changes  in  the  polymer  content 
are  also  obtained.  Dividing  Equation  8  by  Equation  9  gives 
the  expression 

4i  dpi 


Equation  10  divided  by  Equation  9  results  in 


Jfppl 
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Equation  11  divided  by  Equation  9  and  Equation  12  divided  by 
Equation  9  give 

iM  _  rSitP](Uj/r5lvlS]/rJirlPj^ 
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Equations  20  and  21  combine  to  give 
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mass  balance  gives  (dropping  brackets  for  concentration  in  the  remaining  equat 
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where  M 0  is  h  constant  given  by  the  equation 

M*  =  A*  4  B°  +  Pe  GW) 

A 0 ,  E°  and  Pe  representing  initial  quantities  and  dx  is  some 
increment  of  conversion.  The  integral  composition  of  the 
copolymer  is  calculated  stepwise.  Let 

Uj  --  (  dA/J S)' 

*  (dP/Ag). 

Ob  -  dA'/dZ° 

V.  --  <*P‘/d  6* 

and  P  *  r  q 
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where  x  =  (j  +  l)dx.  The  integral  compositions  at  the 
conversion  x  are 

a  =  X  (  J&N ; 
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T  i  u  s  >  i  h  ]  s  k  i  1 1 t- 1  i  c  mode-  1  ,  \olue‘  for  the-  rate  and 
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to  1  in  this  analysis.  The  reactivity  of  the  pendant  vinyl 
rr,a\  be  different  due  to  influences  from  the  polymer  chain 
and  steric  hindrance  and  has  not  been  determined  prior  to 
analysis  with  this  model.  The  constant  kPP  is  equal  to  the 
y-intercept  of  the  pendant  vinyl  conversion  vs.  monomer 
conversion  curve  for  each  system.  The  slope  of  each  of 
these  curves  is  a  result  of  crosslinking  and  secodary 
cyclization  in  combination  with  the  reactivity  of  the 
pendant  vinyl.  If  it  is  assumed  that 
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Figure  1.  Monomer  conversion  vs.  time  of  four  systems 
investigated  . 

Figure  2.  'H  NMR  spectrum  of  sample  from  chemical  system 
2N  . 

Figure  3.  Enlargement  of  4.0  to  6.5  ppm  region  of  2N  ‘H 
NMR. 

Figure  4.  4.0  to  6.5  ppm  region  of  2N  *H  NMR  spectrum  with 

residual  monomer. 

Figure  5.  Pendant  conversion  vs.  monomer  conversion  for 
systems  IN  (♦),  2N  (O),  and  3N  (V). 

Figure  6.  ‘  Pendant  conversion  vs.  monomer  conversion  for 

systems  2N  (♦  )  and  2T  ( )  . 

Figure  7.  Pendant  conversion  vs.  monomer  conversion  for 
systems  2N  (+  )  ,  2N-50  (<>>  and  2N-25  (■). 

Figure  8.  Process  of  crosslink  and  secondary  cycle 
f  orma  t i on . 
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Introduction 

With  the  increasing  use  of  epoxies,  particularly  in  the  aircraft  and  automotive  industries, 
the  establishment  of  direct  correlations  between  network  structure  and  rheological  properties 
is  important.  Information  concerning  structure-rheology  relations  would  assist  in  the 
selection  of  optimum  parameters  for  reactive  processing  operations.  In  addition,  knowledge 
of  these  relationships  can  provide  a  general  understanding  of  the  polymerization  and  network 
formation  processes. 

The  overall  strategy  for  determining  and  applying  structure-property  relations  for  epoxy 
networks  is  outlined  in  Figure  1.  In  order  to  develop  structure-property  relations,  it  is 
essential  to  know  the  way  the  network  links  together.  This  requires  characterization  of  both 
the  reactants  and  the  linking  reaction.  Reaction  kinetics  are  utilized  in  branching  theory  in 
order  to  predict  average  structural  parameters  before  and  after  the  gel  point.  Once  the 
structure  of  the  network  is  ascertained,  it  will  be  correlated  with  rheological  properties. 
Knowledge  of  these  structure-property  relations  can  then  be  used  to  form  a  basis  for  reactive 
process  models. 

This  paper  discusses  the  use  of  branching  theory  to  determine  the  structural  changes 
which  occur  during  the  cure  of  an  epoxy  resin  with  a  primary  amine.  Expenmental  results 
are  presented  to  test  the  ability  of  these  theoretical  models  to  accurately  predict  certain 
network  properties  at  various  stages  of  curing.  Then  these  structural  changes  are  correlated 
with  the  viscosity  rise  during  the  cure. 


aPresent  Address:  Dept,  of  Electrical  Engineering  and  Computer  Science,  Room 
39-561,  Massachusetts  Institute  of  Technology.  Cambridge.  MA  02139  U.S.A. 


Chemical  System  and  Kinetics 


In  order  to  derive  structure-property  relations  tor  any  polymer  system,  specific 
information  is  needed  concerning  the  composition  ot  the  monomeric  units  and  the  chemical 
reactions  which  link  the  units  together.  However,  commercial  epoxy  resin  formulations  are 
generally  quite  complex.  An  accurate  assessment  of  the  initial  composition  of  the  commercial 
formulation  and  its  cure  kinetics  is  virtually  impossible.  Therefore,  to  begin  to  understand 
structure-property  relations  for  commercial  resins,  it  is  necessary  to  first  investigate  and 
understand  network  formation  in  simpler  model  systems. 

The  model  system  chosen  for  this  study  consists  of  diglycidyl  ether  of  bisphenol  A 
(DGEBA)  cured  with  diaminodiphenyl  sulfone  (DDS).  The  chemical  structure  of  the  epoxy 
and  amine  monomer  are  illustrated  in  Figure  2.  Commercial  grade  DGEBA  was  obtained 
under  the  trade  name  DER  332  (Dow)  and  was  purified  by  recrystallization.  The  epoxide 
equivalent  weight  determined  by  chemical  titration  was  173  g/eq.  The  curing  agent  DDS  is 
available  at  99%  purity  from  the  Aldrich  Chemical  Company. 

In  order  to  use  branching  theory  to  predict  network  structure,  one  needs  specific 
information  concerning  the  cure  mechanism.  It  is  necessary  to  know  which  reactions  occur 
during  the  cure,  and  the  ratio  of  rate  constants  for  these  various  reactions.  Shechter  et  al.(1) 
have  suggested  the  sequence  of  reaction  shown  in  Figure  3  for  an  epoxy  cured  with  a 
primary  amine. 

One  factor  which  must  be  considered  when  investigating  the  formation  of  epoxide-amine 
network  polymers  is  the  difference  between  the  reactivity  of  functional  groups  associated 
with  monomers  which  have  been  reacted  at  one  or  more  sites  and  those  attached  to 
completely  unreacted  monomers.  Once  one  of  the  amino  or  epoxide  groups  has  combined 
with  a  bulky  group,  steric  hindrances  and  mobility  restrictions  may  lower  the  rate  at  which 
the  second  functional  group  reacts.  This  difference  in  reactivity  can  have  a  strong  effect  on 
the  structural  development  during  curing. 

Several  workers(2'6>  have  found  that  the  reactivities  of  the  epoxide  groups  in  DGEBA  are 
independent.  In  addition,  the  reactivity  of  amino  groups  in  diamines  also  appears  to  be 
independent. (7)  However,  substitution  among  the  amine  hydrogens  due  to  steric  effect  is  to 
be  expected;  that  is  the  reactivities  of  the  hydrogens  in  the  primary  (rate  constant  fcj )  and  the 
secondary  (rate  constant  k -,)  amino  groups  may  be  different.  The  existing  results  concerning 
the  actual  magnitude  of  the  substitution  effect  in  aromatic  amines  is  somewhat  controversial. 


The  reported  values  for  DDS  show  that  the  ratio  of  rate  constants  U'j/Aa)  range  from 
2  - .  6 .  7  ■ 9  ^ 


The  third  reaction  (etherification)  is  thought  to  occur  by  ternary  amine  or  base  catalyzed 
mechanisms.  Formation  of  polyether  oligomer  has  not  been  observed  in  the  DGEBA-DDS 
systems  when  the  amine  is  present  in  excess  or  at  stoichiometric  concentrations/ :'I0) 
However,  when  epoxide  is  present  in  excess,  the  hydroxyl  groups  may  become  serious 
contenders  with  the  amines  for  the  consumption  of  epoxy  groups/10,111 

All  three  of  the  reactions  which  occur  in  epoxy-amine  systems  can  be  accelerated  by  the 
presence  of  alcohol,  water  or  other  hydrogen  donors/11  Since  a  hydroxyl  group  is  a  product 
of  the  epoxide-amine  addition,  the  course  of  reaction  will  be  autocatalytic. 

Using  this  information  concerning  the  cure  mechanism,  a  kinetic  model  for  the  formation 
of  various  types  of  units  during  the  polymerization.  The  assumptions  made  in  the  kinetic 
scheme  are  listed  below: 

1 )  The  reactivities  of  the  epoxide  groups  are  independc.  ,:t 

2)  The  reactivities  of  the  amino  groups  are  independent 

3)  The  reactivities  of  the  amino  hydrogens  are  dependent 

4)  No  etherification  occurs 

5)  The  presence  of  hydroxyl  groups  catalyzes  the  arrune  -  epoxy  reaction 

6)  The  ratio  of  rate  constants  does  not  depend  on  conversion 

Suppose  the  reaction  begins  with  A4  moles  of  amine  monomer  and  E->  moles  of 
DGEBA.  Then  the  stoichiometric  ratio  is  defined  as 

r  =  2E2/4A4  (1) 

If  at  time  t  a  fraction  pA  of  the  amino-hvdrogens  has  reacted  and  a  fraction  pE  of  epoxide 
groups  has  reacted,  then 

Pa  =  fPe  (2) 


Furthermore,  assume  at  time  t  that  there  are  A4  0  moles  of  A4  units  with  0  reacted  sites,  A4  j 
moles  of  amine  units  with  1  reacted  site,  .  .  .  ,  A44  moles  with  4  reacted  sites.  Thus,  the 
proportion  of  amine  units  which  have  i  reacted  sites  will  be  defined  as  follows: 


Then  the  overall  extent  of  reaction  of  the  amino-hydrogen  groups  can  be  written: 


4 

PA=  X(i/4)pAl  (41 

i=0 


Using  the  above  assumptions  and  notation,  the  cure  of  DGEBA  with  DDS  can  be 
described  by  a  system  of  differential  equations  (Figure  4).  If  it  is  assumed  that  the  ratio  of 
uncatalyzed  primary  to  secondary  amine  reaction  is  the  same  as  the  ratio  of  the  catalyzed, (3) 
the  system  of  equations  have  an  analytical  solution: 


Pa,i  = 

2x(pA.oy  -  Pa.o> 

(5) 

Pa.2*  = 

x2(-2  PA.oy  +  Pa.o+  PA,oa/2) 

(6) 

Pa, 2'  = 

-2x  PA,oy+ax  Pa.o+2  Pa,o1/2 

(7) 

PA.2  = 

Pa, 2*  +  Pa, 2' 

(8) 

PA.3  = 

x2((a+2)  pA/-a  pA  ()-  (2-oupA01--2  pA  *2+(2-a)  pA^4} 

(9) 

II 

-t 

< 

O. 

1  ‘  Pa.o  '  Pa.i  '  Pa. 2  "  Pa.? 

(10) 

where 

cx  —  k- —  kjk^ 

x  =  1/(1 -a/2) 

y  =  ( 1  +  a/2  )/2 

Through  the  system  of  equations  (51  -  (10),  the  structural  configuration  of  all  the 
monomer  units  dunng  the  cure  is  known  as  a  function  of  the  fraction  of  amine  monomeric 
units  which  are  completely  unreacted  fpA  n),  the  initial  stoichiometric  ratio  of  epoxide  groups 
to  amino-hydrogen  groups  (ri.  and  the  ratio  of  rate  constants  (a).  The  reaction  vector  [pA  j, 
. . . .  pA  4,  pE]  is  used  in  the  probability  model  described  in  the  next  section  for  the  prediction 
of  structural  parameters  dunng  poUmerization. 


Branching  Theory 


As  an  epoxy  cures  with  an  amine,  complex  branched  molecules  are  formed.  The 
development  of  these  three  dimensional  structures  determines  processing  properties, 
including  viscosity,  gel  point,  crosslinking  and  modulus.  Following  the  structural  build-up 
of  the  network  during  polymerization  is  difficult.  Before  gelation,  branched  systems  are 
characterized  by  a  wide  molecular  weight  distribution.  Experimental  determination  of  this 
distribution  is  difficult.  After  the  gel  point  pan  of  the  reacting  system  is  incorporated  into  an 
insoluble  gel  with  an  infinite  weight  average  molecular  weight.  Very  few  methods  exist 
which  can  characterize  the  gel  structure.  Clearly  theoretical  relations  capable  of  predicting 
network  structure  given  the  starting  materials  and  the  extent  of  reaction  would  be  valuable. 

Several  statistical  methods  have  been  developed  to  calculate  relations  between  the  extent 
of  reaction  and  molecular  parameters.  Flory(12)  and  Stockmayer(13)  established  the  basic 
theoretical  concepts  of  a  tree-like  model  to  describe  the  branching  of  a  polymer  during 
curing.  However,  this  technique  proved  only  useful  for  very  simplistic,  ideal  reaction 
mechanisms.  More  recently,  two  other  approaches  have  been  developed  to  extend  the  basic 
concepts  set  forth  by  Flory  and  Stockmayer.  The  first  approach  by  Gordon  and 
co-workers(14,15)  uses  the  theory  of  stochastic  branching  processes  and  requires  deriving 
probability  generating  functions.  Macosko  and  Miller1'6'17'  have  developed  an  alternative 
approach  for  predicting  network  structure  which  is  based  on  the  recursive  nature  of  the 
branching  process  and  elementary  taws. 

Using  the  recursive  approach  and  the  assumptions  listed  in  proceeding  section 
concerning  the  reaction  mechanism,  an  expression  has  been  derived(17)  for  the  weight 

average  molecular  weight  as  a  function  of  the  reaction  vector  [pA  j . Pa  4’  PeI-  and  the 

stoichiometric  ratio  of  epoxy  to  amino-hvdrogen  groups  (r).  Figure  5  shows  the  predicted 
change  of  the  weight  average  molecular  weight  with  increasing  conversion  of  epoxy  groups 
for  three  different  stoichiometric  ratios.  The  reaction  vector  is  dependent  on  the  ratio  of  rate 
constants  of  the  secondary  amine  to  primary  amine  reaction,  and  in  Figure  5  a  is  assumed  to 
be  .2. 

The  gel  point  is  defined  as  the  conversion  at  which  the  weight  average  molecular  weight 
diverges  and  an  infinite  network  begins  to  form.  As  is  illustrated  in  Figure  5,  the  molecular 
weight  builds  very  slowly  until  near  the  gel  point.  The  conversion  at  gelation  is  highly 
dependent  on  the  initial  ratio  of  epoxy  groups  to  amino- hydrogen  groups.  An  expression  can 
be  derived  for  the  gel  conversion  of  epoxy  groups  (Figure  6),  by  setting  the  molecular 
weight  equal  to  infinity  and  solving  for  the  conversion  as  a  function  of  the  stoichiometric 


ratio.(17)  There  is  a  point  at  which  the  ratio  of  epoxide  to  amino-hydrogen  groups  is  so  small 
that  the  system  will  not  gel.  This  is  known  as  the  critical  stoichiometric  ratio.  Under  the 
condition  that  a  =  .2,  the  critical  stoichiometric  ratio  is  predicted  to  be  .42. 

The  recursive  approach  can  also  be  used  to  predict  post-gel  properties.  Using  the 
preceding  assumptions  concerning  the  epoxy-amine  cure  reaction,  an  expression  has  been 
derived  for  predicting  the  weight  fraction  of  soluble  material.  ■  Figure  7  shows  the 
predicted  change  of  sol  fraction  as  a  function  of  epoxide  group  conversion,  assuming  a  =  .2. 

Another  parameter  which  can  be  derived  from  the  recurs i%  e  theory  is  the  weight  average 
molecular  weight  of  the  longest  chain  which  passes  through  a  randomly  chosen  point  on  the 
molecule.  This  structural  parameter  reflects  the  hydrodynamic  size  of  a  branched  polymer 
and  thus  has  been  found  useful  in  correlating  the  viscosity  rise  during  polymerization. (l81 
Miller  et  al/19)  have  derived  the  longest  linear  chain  for  a  chemical  system  which  has  equal 
reactivity  of  functional  groups  of  the  same  type,  and  w here  the  monomenc  molecular  weight 
of  the  crosslinker  is  negligible  compared  to  the  weight  of  the  other  monomer.  This 
derivation  has  been  extended  to  include  the  substitution  effect  and  the  monomer  weight  of  the 
amine  crosslinker/20)  Figure  8  compares  the  rise  of  the  weight  average  molecular  weight 
with  the  molecular  weight  of  the  longest  linear  chain  at  balanced  stoichiometry. 

Figures  5  -  8  assume  that  the  ratio  of  rate  constants  for  the  secondary  amine  to  primary 
amine  reaction  is  .2.  However,  as  discussed  in  the  previous  section,  the  magnitude  of  the 
substitution  effect  in  the  DDS  system  is  controversial.  The  literature  reports  values  ranging 
from  .2  -  .6.  In  systems  with  complicated  cure  chemistry,  it  is  frequently  difficult  to 
ascertain  specific  values  for  the  rate  constant  ratios  of  the  various  reactions  which  occur 
during  polymerization.  Hence,  it  is  important  to  consider  the  sensitivity  of  the  predictions 
for  structural  parameters  to  variations  in  rate  constant  ratios.  Figure  9  examines  the  change 
in  gel  conversion  with  the  magnitude  of  the  substitution  effect.  At  balanced  stoichiometry, 
when  a  =  .2  the  predicted  gel  conversion  is  .61,  while  the  gel  conversion  is  .59  when  a  = 
.6. 

Experimental 

The  goal  of  this  study  is  to  develop  relationships  between  rheological  properties  and  the 
build  up  of  network  structure  during  polymerization.  Therefore,  it  is  first  necessary  to  test 
the  ability  of  the  branching  expressions  to  predict  molecular  parameters.  Then  these 
expressions  can  be  combined  with  rheological  and  kinetic  data  to  form  structure-property 
models. 
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As  described  in  previous  section,  the  recursive  method  is  capable  of  predicting 
conversion  at  gelation  as  a  function  of  the  stoichiometric  ratio  of  epoxide  to  amino- hydrogen 
groups.  This  relationship  can  also  be  empirically  determined.  Samples  of  given 
stoichiometric  ratios  were  cured  in  narrow  test  tubes  immersed  in  a  177:C  oil  bath.  At 
various  reaction  times,  the  samples  were  removed  from  the  oil  bath  and  rapidly  plunged  into 
an  ice  bath  to  stop  the  reaction.  The  samples  were  then  removed  from  the  test  tubes  and 
extracted  with  tetrahydrofuran.  The  gel  time  was  designated  as  the  cure  time  at  which  the 
sample  just  can  not  dissolve  in  the  solvent.  The  conversion  of  epoxide  groups  at  this  time  is 
determined  by  differential  scanning  calorimetry  (DSO.l20)  These  results  are  plotted  along 
with  the  predictions  of  the  recursive  theory  (assuming  cc=.  2)  in  Figure  10.  Good  agreement 
is  seen  between  the  experimental  points  and  the  recursive  theory. 


The  critical  stoichiometric  ratio  predicted  by  the  branching  theory  can  be  verified  by 
curing  samples  with  different  stoichiometric  ratios  and  performing  extraction  experiments  to 
determine  gel  content.  Several  DGEBA-DDS  samples  were  prepared  with  the  stoichiometric 
ratio  of  epoxide  to  amino-hydrogen  groups  ranging  from  .2  to  .6.  The  samples  were  placed 
under  nitrogen  in  a  177°C  oven  for  approximately  24  hours,  and  then  were  post-cured  at 
230°C  for  8  hours.  The  samples  were  then  crushed,  placed  in  cellulose  thimbles,  and 
extracted  for  a  week  in  tetrahydrofuran.  The  critical  stoichiometric  ratio  is  the  largest  ratio  of 
epoxide  to  amino-hvdrogen  groups  where  the  polymer  is  still  completely  soluble.  The 
critical  stoichiometric  ratio  for  DGEBA-DDS  was  determined  to  be  .40.  Comparison 
between  this  value  and  the  predicted  value  from  the  recursive  theory  is  shown  in  Figure  10. 


The  sol  fractions  of  samples  cured  beyond  the  gel  point  were  also  measured  for  the 
DGEBA  system.  Samples  were  cured  at  various  time  in  a  177°C  oil  bath.  The  reaction  was 
stopped  by  quenching  in  an  ice  bath.  Pre-weighed  samples  were  crushed,  placed  in 
extraction  thimbles,  and  extracted  in  tetrahydrofuran  for  several  w-eeks,  with  the  solvent 
exchanged  periodically.  Then  the  samples  were  dried  under  vacuum  until  no  further  weight 
loss  could  be  detected.  The  final  sample  weight  was  obtained  and  the  weight  fraction  of 
solubles  was  calculated.  The  conversion  of  epoxide  groups  at  a  given  cure  time  was 
determined  by  DSC.  Figure  1 1  shows  excellent  agreement  between  theoretical  and 
experimental  values  obtained  for  DGEBA-DDS  at  balanced  stoichiometry. 

Correlations  between  structure  and  rheological  properties  can  be  established  by 
simultaneously  measuring  viscosity  and  conversion  during  cure.  The  steady  shear  viscosity 
nse  for  the  DGEBA-DDS  system  has  been  monitored  using  the  parallel  plate  geometry'  on  the 
Rheometncs  System  Four  Rheometer  The  cure  was  conducted  at  177CC  under  a  nitrogen 
environment.  At  specific  times  during  the  cure,  the  reaction  was  stopped  in  the  rheometer  by 


quenching  in  liquid  nitrogen.  Sample  was  then  removed  from  the  rheometer  and  analyzed 
for  conversion  by  DSC.  Thus,  the  viscosity  rise  can  be  obtained  as  a  function  of  conversion 
(Figure  12).  Using  the  recursive  theory,  the  molecular  weight  of  the  longest  linear  chain  can 
be  determined  at  any  conversion.  Thus  in  Figure  13  the  viscosity  is  plotted  as  a  function  of 
molecular  weight  of  the  longest  linear  chain.  As  seen,  the  viscosity  data  for  different 
stoichiometric  ratios  collapses  onto  a  single  curve. 

It  has  been  shown  that  the  viscosity  for  many  undilute  polymers  with  chain  lengths  less 
than  300  -  500  main  chain  atoms  is  dependent  on  the  friction  factor  g  and  molecular  weight 
MW; 

na!;M  (11) 

The  friction  factor  is  governed  by  local  intrachain  and  interchain  forces  between  neighboring 
segments  in  the  polymer  melt.  Its  magnitude  is  dependent  on  the  particular  polymer  system. 
A  general  expression  has  been  derived  for  the  friction  factor  as  a  function  of  material 
properties  of  the  polymer  system(21): 

£  =  C0  exp[B/(fg  +  a,  (T  -  Tg))]  (12) 

where  cq  is  the  liquid  expansion  factor,  f  is  the  fractional  void  or  free  volume  at  the  glass 
temperature  Tg,  B  is  related  to  the  fractional  void  volume  required  to  be  in  the  vicinity  of  a 
segment  for  that  segment  to  make  a  jump,  and  Q)"1  is  an  inherent  jump  factor. 

In  order  to  correlate  the  viscosity  rise  with  molecular  weight,  it  is  necessary  to  calculate 
the  friction  factor  at  any  point  in  the  polymerization.  Through  the  examination  of  several 
polymer  systems.  Berry  and  Fox(21)  found  that  the  inherent  jump  factor  and  the  ratio  B/fg  are 
generally  constant  for  any  system.  However,  Tg  and  to  a  lesser  extent  ccj,  depend  on 
molecular  weight  until  somewhat  longer  chain  lengths  (i.e.  greater  than  ca.  500)  are  attained; 
then  they  too  become  essentially  molecular  weight  invariant.  The  ratio  10‘3B/  a,  is  in  the 
interval  1.5  to  3.0  with  B  increasing  with  the  size  of  any  side  groups.  In  order  to  calculate 
the  friction  factor  for  the  DGEBA-DDS  system  it  will  be  assumed  that  the  ratio  10'3B/  a,  is  a 
constant,  equal  to  2.25.  The  glass  transition  temperature  is  more  sensitive  than  this  ratio  to 
changes  in  molecular  weight  during  polymerization.  According  to  the  approach  of  Adabbo 
and  Williams(22)  using  the  DiBenedetto123’  equation,  the  glass  transition  temperature  increases 
with  increasing  extent  of  conversion  as  follows: 


where  Ex/Em  is  the  ratio  of  lattice  energies  for  crosslinked  and  uncrosslinked  polymer, 
Fx/Fm  is  the  corresponding  ratio  of  segmental  mobilities,  Xg  is  the  conversion  at  the  glass 
transition  temperature  and  Tg0  is  the  glass  transition  temperature  at  no  reaction.  The  ratios 
Ex/Em  and  Fx/Fm  have  been  determined  by  fitting  equation  13  to  conversion  versus  glass 
transition  temperature  data.  For  the  DGEBA-DDS  system  at  balanced  stoichiometry,  the 
following  constants  have  been  obtained  for  the  DiBenedetto  equation-24); 

Ex/Em=  1-21 
=  0.74 


Using  the  typical  material  constants  obtained  by  Berry  and  Fox  (1968)  and  the 
DiBenedetto  relationship  for  glass  transition  temperature,  the  friction  factor  £  can  be 
estimated  for  balanced  stoichiometry  as  a  function  of  conversion  of  epoxide  groups.  In 
Figure  14,  viscosity  is  plotted  against  the  product  of  the  estimated  friction  factor  and  the 
molecular  weight  of  the  longest  linear  chain  for  the  DGEBA-DDS  system  at  balanced 
stoichiometry.  Linear  regression  of  the  data  gives  a  slope  of  1.08  with  a  regression 
coefficient  of  0.98.  Thus,  the  viscosity  rise  for  the  DGEBA-DDS  system  at  balanced 
stoichiometry  fulfills  the  relation  of  being  dependent  on  a  friction  factor  and  molecular  weight 
to  the  first  power.  A  possible  explanation  for  the  curvature  shown  in  the  data  points  in 
Figure  14  may  be  a  result  of  the  assumption  that  B/  dj  is  a  constant  during  polymerization. 
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Figure  3.  Reaction  mechanism  for  an  epoxy  cured  with  a  primary  amine 
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Figure  4.  Kinetic  scheme  for  the  DGEBA-DDS  system. 


MOLECULAR 

WEIGHT 

10- 


10’ 

0  0  0.2  0.4  0.6  0.8  l.G 

CONVERSION 

Figure  5.  Recursive  model  predictions  for  molecular  weight  as  a  function  of  conversion 
at  three  different  stoichiometric  ratios.  For  the  model  predictions,  the  ratio  of 
the  secondary  amine  to  primary  amine  (a  )  is  assumed  to  be  .2. 


Figure  6.  Recursive  model  prediction  for  the  gel  conversion  of  epoxide  groups  as 
a  function  of  the  initial  stoichiometric  ratio  of  epoxide  groups  to  amino- 
hydrogens.  For  the  model  predictions,  the  ratio  of  the  primary  to 
secondary  amine  reaction  (a )  is  assumed  to  be  .2. 
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Figure  7.  Recursive  model  predictions  for  the  weight  fraction  of  soluble  material 
versus  the  conversion  of  epoxide  groups  at  three  stoichiometric  ratios. 
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Figure  8.  Comparison  of  the  weight  average  molecular  weight  and  the  molecular  weight 
of  the  longest  linear  chain  at  gelation. 
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Figure  9.  Conversion  of  epoxide  groups  at  gelation  versus  the  rate  constant  ratio 
of  the  secondary  amine  to  the  primary  amine  reaction  for  three  different 
stoichiometries. 
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Figure  10.  Gel  conversion  of  epoxide  groups  versus  stoichiometric  ratio.  Predictions 
of  the  recursive  theory  are  indicated  by  the  solid  line.  Experimental  gel 
conversion  values  (  ■  )  were  obtained  by  DSC. 
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Figure  13. 


Viscosity  rise  versus  molecular  weight  of  the  average  longest  linear  chain  at 
stoichiometric  ratios  ranging  from  .6  to  2.0.  The  molecular  weight  of  the 
longest  linear  chain  is  calculated  using  the  recursive  theory  assuming  a  =  .2. 


Figure  14.  Viscosity  rise  versus  the  product  of  the  friction  factor  and  the  molecular  weight 
of  the  average  longest  linear  chain  for  the  DGEBA-DDS  system  at  balanced 
stoichiometry. 


